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A spin-tunnel-junction based on manganites, with Lai-ajSr^MnOs (LSMO) as ferromagnetic 
metallic electrodes and the undoped parent compound LaMnOs (LMO) as insulating barrier, is 
here theoretically discussed using double exchange model Hamiltonians and numerical techniques. 
For an even number of LMO layers, the ground state is shown to have anti-parallel LSMO magnetic 
moments. This highly resistive, but fragile, state is easily destabilized by small magnetic fields, 
which orient the LSMO moments in the direction of the field. The magnetoresistance associated 
with this transition is very large, according to Monte Carlo and Density Matrix Renormalization 
Group studies. The influence of temperature, the case of an odd number of LMO layers, and the 
differences between LMO and SrTiOa as barriers are also addressed. General trends are discussed. 

PACS numbers: 73.90.-|-f, 71.10.-w, 73.40.-c, 73.21.-b 



I. INTRODUCTION 

The study of strongly correlated electronic systems 
(SCES) continues attracting the attention of the Con- 
densed Matter community. These materials present com- 
plex phase diagrams that illustrates the competition 
which exists among phases with very different physical 
properties, such as c?-wave superconductivity, antiferro- 
and ferro-magnetic order, charge- and orbital-order, mul- 
tiferroic behavior, and several others. Moreover, this 
complexity and phase competition lead to self-organized 
nano-scale inhomogeneities which are believed to gener- 
ate giant responses, as in the famous colossal magne- 
toresistance (CMR) effect of the Mn-oxides known as 
manganites.^ 

Recently, a new procedure to study oxide SCES has 
been proposed. It involves the artificial creation of oxide 
multilayers with atomic-scale accuracy at the interface, 
via the use of techniques such as pulsed-laser deposition^ 
Potentially, these structures can have properties very dif- 
ferent from those of the building blocks. One of the rea- 
sons for this expectation is that a transfer of charge could 
occur between the constituents leading, for example, to 
the stabilization of a metal at the interface between two 
insulators.^ The creation of novel two-dimensional states, 
as well as the possible applications of oxide multilayers 
in the growing field of oxide electronics, has given con- 
siderable momentum to these investigations. 



A. Spin-Tunnel-Junctions 

The development of the above mentioned accurate ex- 
perimental techniques for the construction of oxide mul- 



tilayers with atomic precision can have implications in 
the study of spin-tunnel-junctions, ■^■'^ introducing a bet- 
ter control of their properties. These structures con- 
sist of two ferromagnetic (FM) metallic electrodes, sep- 
arated by a thin insulating barrier. The resistance of 
this device depends on the relative orientation of the 
electrodes' magnetizations. The tunneling magnetoresis- 
tance (TMR) is usually defined via the difference in re- 
sistances between the anti-parallel and parallel arrange- 
ments of the electrodes' magnetic moments. Half-metals, 
such as Lai_a;Sra;Mn03 (LSMO), with an intrinsic nearly 
full magnetization, are ideal for these devices*^ 

In the context of magnetic tunnel junctions, very in- 
teresting results were reported by Bowen et al,^ using a 
LSMO/STO/LSMO trilayer. The Sr concentration was 
1/3 for Lai_a;Sra;Mn03 (LSMO), and STO represents the 
insulator SrTiOa. A huge TMR ratio of more than 1800% 
was observed at very low temperatures 4 K, showing the 
advantages of using half-metallic LSMO as a ferromag- 
netic electrode in the junctions. However, in the same 
investigations it was reported that the large TMR sur- 
vived only up to ^270 K, lower than the Curie temper- 
ature of LSMO (a;= 1/3), which is ~370 K. It was argued 
that the deterioration of the ferromagnetism near the 
LSMO/STO interface ("dead layer") could be causing 
this TMR reduction. Later, Yamada et alX addressed 
this problem by comparing the STO/LSMO interface 
with others, such as LAO/LSMO or STO/LMO/LSMO, 
where LAO stands for LaAlOa and LMO for LaMnOa. 
Those authors found that the magnetic behavior of 
LAO/LSMO and STO/LMO/LSMO are much better 
than STO/LSMO in the sense that no dead-layer was 
found, opening a new path toward LSMO-based TMR 
junctions operating at room temperature. 
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B. Proposed Main Idea 



In this paper, an alternative setup is proposed for a 
manganite trilayer system which is expected to have a 
very large magnetoresistance (MR) at low temperatures, 
at least according to modeling calculations reported be- 
low. The proposed geometry, and main idea behind its 
performance, is presented in Fig. [TJ The new system is 
made out entirely of manganite materials, with different 
hole doping concentrations.'^ Using all manganites may 
help in the interfacial contact between the components 
due to their similar lattice spacings. More specifically, 
in Fig. [1] (a) a trilayer system is represented. It contains 
two hole-doped Mn-oxides, such as LSMO, and a central 
region made out of the undoped parent compound LMO. 
It is well known that LSMO is ferromagnetic and half- 
metallic for sufficiently large hole doping, while LMO is 
an A-type antiferromagnetic (AF) insulator^ The arrows 
in Fig. [T] (a) represent schematically the expected spin 
orientations in a one-dimensional arrangement, for sim- 
plicity. The most interesting detail of Fig. [T] (a) is the rel- 
ative orientation of the spins between the metallic leads. 
For an even number of layers in the central LMO region, 
the ferromagnetic moments of the leads are anti-parallel. 
This is expected to cause a very large resistance at low 
temperatures, since the carriers moving from one lead to 
the other not only must tunnel through the central insu- 
lating barrier, but in addition the spin species which can 
travel in one lead is blocked by the other. However, the 
anti-parallel configuration, which is mediated by the cen- 
tral region, is not strongly pinned in this arrangement: 
it is only a weak antiferromagnetic effective interaction 
which produces the ground state with anti-parallel leads' 
moments. Thus, relatively small magnetic fields can ren- 
der the ferromagnetic moments parallel, substantially re- 
ducing the resistance. All these intuitive ideas will be 
substantiated via model calculations, described below. 
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FIG. 1: (Color online) Schematic representation of a 
one- dimensional (ID) double exchange (DE) model for a 
LSMO/LMO/LSMO trilayer with even (a) and odd (b) num- 
bers of LMO sites. Here, LMO and LSMO represent regions 
where the positive background charge densities correspond 
to an antiferromagnetic (AFM) insulator and ferromagnetic 
(FM) metal, respectively, when in bulk form. Circles stand 
for lattice sites where electrons can move, and arrows in- 
dicate localized spins. The trilayer with even (odd) num- 
ber of LMO sites gives rise to anti-parallel (parallel) align- 
ment of the magnetic moments in the left and right ends 
of LSMO. (c) Schematic representation of a ID model for 
a LSMO/STO/LSMO trilayer. Here, LSMO is described by 
the DE model as before, while STO is modeled by a tight- 
banding model with no localized spins. Vm introduced in the 
central STO region is a band off-set site-potential mimicking 
the work function difference between LSMO and STO. More 
details can be found in the text. 



C. Main approximations 



For an odd number of layers in the central LMO region, 
the expected spin arrangement is shown in Fig. [T] (b) . In 
this case the magnetic moments of the leads will be paral- 
lel to one another, and the resistance will not be as large 
as for the configuration shown in Fig. [1] (a). However, 
the LMO region still provided a tunneling barrier, and 
the performance of this case will be shown below to be 
quite acceptable, at least within modeling calculations. 

Finally, Fig. [TJc) contains a crude representation of a 
more standard trilayer device.- While the leads are still 
representing LSMO, the central region is now a band in- 
sulator, as it occurs in the much employed case of STO. 
The main focus of our effort will be the cases shown in 
Fig. [T] (a) and (b) , but some results for the setup shown 
in Fig. [T] (c) will also be discussed for completeness, and 
to clarify qualitative trends. 



Before proceeding to the presentation of the results, 
some of the theoretical approximations used must be 
clearly expressed for the benefit of the reader. With this 
paper, our main intention is to motivate experimental 
groups to consider the materials and setups proposed in 
our study, involving a magnetically active manganite bar- 
rier, as opposed to a magnetically-inert band-insulator, 
as the widely used STO. However, it must be clearly 
stated that our calculations are qualitative at best and 
should be considered only as a guidance to understand 
the intuitive picture presented here. 

A variety of effects are not taken into account in this 
investigation (and actually these effects cannot be taken 
into account accurately within the current status of nu- 
merical simulations and Hamiltonian modeling). (1) For 
instance, the lattice and orbital reconstructions are not 
incorporated, but only the electronic reconstruction is 
considered. In other words, the atomic positions are con- 
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sidered to be rigid here. It would be the task of sophis- 
ticated ab-initio simulations to consider how the lattice 
relaxes at the interface between LMO and LSMO, and 
such future calculations are certainly strongly encour- 
aged. (2) A second subject which is not addressed in 
this work is the issue of the infamous "dead layer" at 
the barrier/electrode interface, already briefly described 
in Sec. II Al Dead layers affect the performance of sev- 
eral trilayer devices. We have attempted to mimic this 
dead layer altering by hand the chemical potential in 
the vicinity of the interface, but the results were not 
sufRciently satisfactory to be described in this paper. 
Thus, this problem is left for future efforts. (3) Another 
topic which is only briefly discussed is the influence of 
anisotropics: in most of the simulations below there is no 
"easy" axis introduced for the magnetization orientation. 
While for manganites in bulk form this is a reasonable 
assumption, in thin films it is known that in-plane mag- 
netic anisotropics are induced by epitaxial strain.^ These 
anisotropics are the reason behind the abrupt changes 
in resistances observed when varying magnetic fields in 
magnetic tunnel junctions (see, for instance. Fig. 2 in 
Ref. H). While in our numerical studies reported be- 
low the simulations are carried out mainly with isotropic 
Heisenberg spins, a more proper analysis would have 
needed anisotropic terms, rendering these spins Ising like. 
(VJ A related topic involves the experimentally observed 
differences between the upper and lower electrodes in a 
trilayer junction, since the area of the bottom electrode 
is typically much wider than that of the upper onei^^ 
This stabilizes anti-parallel orientations of the electrodes' 
magnetic moments at some magnetic fields. In our sim- 
ulations, both electrodes are perfectly equivalent. How- 
ever, for the reasons already mentioned in the previous 
subsection, we do stabilize a similar anti-parallel mag- 
netic moment arrangement via the use of LMO as bar- 
rier, with an even number of sites and in zero external 
magnetic field. (5) Finally, a practical assumption in our 
investigations is the focus on a one-dimensional (ID) spin 
arrangement, described by the realistic double-exchange 
model for manganites (albeit restricted to just one or- 
bital). The restriction to a ID configuration is needed 
for the numerical studies to be accurate; higher dimen- 
sional arrangements would have increased so much the 
CPU time that a careful analysis would have been im- 
possible. Note that contrary to manganite bulk studieSfi 
where the analysis of small two-dimensional clusters is 
possible, here we will be carrying out an iterative loop 
to solve Poisson's equation which regulates the charge 
transfer between materials (see SecHIl. At each step of 
the iterative process, an entire Monte Carlo (MC) pro- 
cess or Density Matrix Renormalization Group (DMRG) 
sweep is carried out (see details below), thus increasing 
substantially the computer requirements. Then, the re- 
striction to a ID geometry is caused by the CPU re- 
sources available. However, in the description of results 
below, we have focus on qualitative aspects which are ex- 
pected to be robust, and we strongly believe that they 



will survive the increase in dimensionality. For example, 
the large magnetoresistance at low temperatures of the 
LSMO/LMO/LSMO setup, with an even number of lay- 
ers for LMO, is believed to occur in any dimension of 
interest. 



D. Organization 

The organization of the paper is the following. In 
Sec. ini the model Hamiltonians and numerical meth- 
ods are described. The focus is on Monte Carlo and 
DMRG techniques. The main results are presented in 
Sec. Illli which correspond to the arrangements schemat- 
ically described in Fig. [T] (a) and (b), showing that the 
magnetoresistance is large in these setups. Both classical 
and quantum localized t2g spins are used. In Sec. llVi 
the influence of temperature is analyzed. As observed in 
some experiments,^ it is found that the large MR effect 
quickly deteriorates with increasing temperature. The in- 
fluence of anisotropics is also studied as a possible cure to 
this problem. In Sec. |Vl for completeness, results of the 
modeling of a trilayer involving a band insulator (such 
as STO) as the barrier, instead of LMO, are reported. 
Conclusions are given in Sec. IVII 

II. MODEL AND METHODS 

A. Model Hamiltonian and MC methods 

To crudely model LSMO/LMO/LSMO trilayers 
(Fig. [T]), the one-orbital double exchange (DE) model on 
a ID lattice will be used: 

H = -t5]5](4,c,,, + H.c.) +^(/)(^)n, 

-Jh^s, -S, -^hext-M,, (1) 

i i 

where c| ^ is the creation operator of an electron at 
site i with spin (t{—'\,1). The first summation of 
runs over nearest neighbor pairs of sites i and j. The 
number operator is rii — ^^c\ ^Ci,^^ cj){i) is the elec- 
trostatic potential (discussed below), and the spin op- 
erator of the electron is — X^a./J 4,q i.^)ai3 ^i,/? (here 

= (fa;, <JyT<Jz)'- Pauli matrices). Si is the classical spin, 
widely used to represent the localized i2g spins, with 
|Si| = l. Jh is the Hund's rule coupling, hcxt is an ex- 
ternal magnetic field, and Mi=i(si + 38^) is the total 
magnetic moment at site i. Hereafter, t is set to be 1 as 
an energy unit. The number of sites in the left, central, 
and right regions of the system is denoted by L''^\ L'-'^\ 
L(^\ respectively. 

To study the electronic properties of hetero-structured 
systems, it is crucial to include the cation ions in the 
model (to consider the charge neutrality condition) and 
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take into account long-range Coulomb interactions be- 
tween electrons and cation ions. In this study, the long- 
range Coulomb interactions are considered within the 
Hartree approximation through Poisson's equationi^i 

V^<p{^) = -c^[{n^)-n+{^)], (2) 

where is the positive background charge mimick- 

ing the cation ions. Here in this paper, n^{i)—n^^'^ 
is set to be uniform within each layer [r=L (left re- 
gion), C (central region), and R (right region)], with 
a value determined by the charge neutrality condition. 
The parameter a=e^ /ea (e: electronic charge, e: dielec- 
tric constant, a: lattice constant) is the strength of the 
Coulomb interactions, considered as a free parameter in 
the model. To solve Eq. ([2]), the symmetric discretiza- 
tion of Poisson's equation is used, which in ID becomes 
V'^ (/){{)— 4i{i + 1) — + — 1), with boundary con- 

ditions (f>{i)~0 for sites outside of the system. 

In the following sections, results for zero temperature 
(T) as well as finite temperatures are reported. For 
the finite temperature calculations, the standard grand- 
canonical Monte Carlo simulation is used^^ with the 
chemical potential /i adjusted such that the total num- 
ber of electrons are equal to the total background positive 
charge n+{i). For the zero temperature calculations, 
a full optimization of {S^} using the Broyden-Fletcher- 
Goldfarb-Shanno method^'^ is performed for a fixed num- 
ber N of electrons, i.e., N = n+li) (canonical ensem- 
ble). We have found that the Monte Carlo simulations 
at very low temperatures produce almost identical results 
as those calculated using the zero temperature canonical 
method, although at a considerably larger cost in CPU 
time. Thus, at T=0 it is advantageous to use the op- 
timization method. Note that due to the need to solve 
Poisson's equation iteratively, a Monte Carlo simulation 
or optimization procedure has to be carried out at each 
step, increasing substantially the computer time as com- 
pared with more standard simulations of bulk systems. 
The calculation of the conductances is carried out us- 
ing the Landauer formalism," as extensively explained 
in previous reports.^'"^^ 

In the real LSMO/STO/LSMO trilayer systems, which 
will also be briefly modeled in this manuscript, it is 
important to notice that there exists band offsets (due 
to work- function differences) between LSMO and STO, 
which are typically of the order of a few eVs,^^ although 
its precise value is difficult to find experimentally.^^ In 
this paper, this band offset is treated as a parameter de- 
scribed simply by the addition of a site potential term 
Vm J2icc Hamiltonian [Fig. [1] (c)]. Moreover, a 

simple tight-binding model without Jh in Eq. ([T]) is used 
to model the STO barrier. Omitting Hubbard-typc in- 
teractions might be justified by the fact that the number 
of electrons in the central region (STO) is very small, 
as shown later. And for LSMO and LMO, neglecting the 
Hubbard term is justified by the large value of the Hund's 
coupling which by itself prevents double occupancy, as 
widely discussed beforcji 



B. Quantum localized spins and DMRG 

Results using quantum spin 1/2 for the localized spins 
Si in Eq. ([T|) are also presented in this paper, and they 
are compared with results of the classical spin simula- 
tions. In real manganites t2g spins are 3/2. However, 
using spin 1/2 much simplifies the computational task 
due to the reduction in the size of the Hilbert spaccA^ 
To study the ground state properties of Eq. ^ with 
Si = 1/2, e.g., the charge density distribution, the stan- 
dard DMRG algorithm,^* embedded in a self-consistent 
iterative procedure to solve Poisson's equation for the 
long-range Coulomb potential, is used. 

More specifically, starting from an initial electrostatic 
potentials (f>{i) {i — 1, . . . , L, here L is number of total 
sites), we first make two sweeps for "warming" before be- 
ginning the self-consistent calculations for the long-range 
Coulomb interactions. Then, for the next 10 sweeps, the 
electrostatic potentials 0(i) are updated by solving Pois- 
son's equation [Eq. ^] at each sweep. In general, the 
convergence of the self-consistent procedure is very slow. 
Therefore, here we perform an extrapolation of for 
each i, using the values calculated during those 10 sweeps, 
to an infinite number of sweeps. By comparing fully self- 
consistent calculations, we have found that generally this 
extrapolation scheme gives a reasonable results. Finally, 
the obtained potential (p{i) is plugged back in Eq. ([T]), 
and 8 more sweeps are performed to calculate the ground 
state. The results reported below are obtained by retain- 
ing Af=350 states, and the truncation error in the worst 
case is of order 10~^. All DMRG calculations are done 
at T = 0, and Jh=36 is used^^^ 

As a check of our DE model code, we have compared 
the case of a total z-projection of the spin Sz = 5'max=26, 
on a lattice of L=32 sites and with a total number of con- 
duction electrons N = 20, with the results of the spin- 
less fermion model which is computed with a previously 
prepared DMRG code for the Hubbard model, written 
completely independently from the present DE code. The 
potential added in the spinless fermion calculation 
is the one extracted from the DE run. The differences in 
energies are found to be of order 10"^. 

To study transport properties, the time-dependent 
DMRG technique is used.— The electrostatic potentials 
(j>{i) are fixed to the ones obtained above, and a small 
bias potential AT^^O.Ol is applied at time t=0, which 
triggers a time-evolution of the systemi^ The bias po- 
tential is applied only on a few sites at the edges of the 
system, and the current as a function of t is measured 
on the two links connecting the central region to the two 
outer portions of the system. The amplitude of J{t), 
defined as the average of these two currents, scales ap- 
proximately linearly with the number of sites on which 
the bias potential AV^ is applied. As it is well known^^ 
J{t) follows an oscillatory evolution with time due to the 
open boundary conditions used in the DMRG process 
(electrons cannot leave the system). Thus, to obtain a 
measure of the conductance, the average of J{t)/ AV over 
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the first half-period of the oscillation is considered i2i 

III. RESULTS AT ZERO TEMPERATURE 

In this section, the main results obtained in our nu- 
merical simulations will be described. 



A. All-manganite trilayer geometry 
using classical t2g spins 

1. Even number of sites in the barrier 

Figure [2] shows the zero-temperature optimization re- 
sults corresponding to a ID trilayer geometry [Fig.[l](a)] 
consisting of 20 sites in each lead, and 4 sites in the mid- 
dle. The positive charge, regulated by n+, is 0.65 at 
each lead and 1.0 in the middle (in e=l units). The 
strength of the long-range Coulomb interaction is chosen 
as a=1.0. Figure dja) shows the converged local elec- 
tronic density n{i) vs. the site location i, along the 44 
sites chain. In the leads, the electronic density closely 
matches the expected result 0.65. The oscillations at the 
end of the chain near i—1 and 44 are due to Friedel os- 
cillations. In the 4-sites center, 2 of the sites have n{i) 
very close to 1, while the other 2 have a smaller density 
due to the charge-transfer effect of the long-range forces. 
The overall electronic density profile is reasonable and in 
agreement with qualitative expectations. 

The lower panel Fig. [^Jb) contains the most impor- 
tant result for this particular geometry. There, the ori- 
entations of the classical t2g spins are shown via the 
spherical coordinates angle Lp, as a function of position i. 
Note first that each lead has classical spins polarized 
in a ferromagnetic state, namely all sites have approx- 
imately the same Lp, as expected from the well-known 
double-exchange mechanism;21 This ferromagnetism at 
the leads is compatible with phase diagrams gathered 
in previous studies However, in the near absence of 
magnetic fields, /icxt= 0.00001, the angle Lp of the two fer- 
romagnetic leads differs by tt, signaling an anti-parallel 
arrangement of magnetic moments between the two ferro- 
magnetic leads (the tiny field was used simply to orient 
the spins along a direction which simplifies the discus- 
sion, and it does not have any other important effect). 
This anti-parallelism is in excellent agreement with the 
expected results based on the introductory discussion: if 
the number of sites in the antiferromagnetic barrier is 
even, then an anti-parallel orientation of ferromagnetic 
moment between the leads should occur. In the barrier 
region, the classical spins are arranged in an antiferro- 
magnetic pattern, since there n(i)~l and an antiferro- 
magnetic spin orientation is preferredj^ 

The results become even more interesting as the mag- 
netic field is increased. In this case, the relative orienta- 
tion of the ferromagnetic moments in the leads changes 
very rapidly, as shown in Fig. [21Ib). For fields as small 
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FIG. 2: (Color online) (a) Local electron density n(i) and 
(b) classical spin orientation Si = {6i,ifii) (0 < 6i < n, 
< LPi < 27r) in spherical coordinates vs. site location i, 
on a trilayer chain with L^^^ = L*^^' = 20 sites in each lead 
and L^'^^ = 4 sites in the center. The ID DE model with 
Jh=8.0 and 0=1.0 is used, and the results are obtained via 
an optimization method at zero temperature. Lead's positive 
background charge is n^'=n^-' =0.65, and a central positive 

(C) 

background charge is =1.0. A magnetic field is applied in 
the direction opposite to a;-direction (x), and the magnitude 
(ftcxt) is shown in the figures. In (b), 6i ~ 7r/2 for all /icxt's. 
As ^ext increases, the spins gradually rotate toward — x in 
the xy plane, and finally they align ferromagnetically. The 
positions of the two interfaces are denoted by vertical dashed 
lines. 



as hcxt— 0.005 (in units of the hopping i), the magnetic 
moments of the ferromagnetic leads are already nearly 
aligned. As shown below, this produces drastic changes 
in the conductance of the ensemble. Results for interme- 
diate values of the magnetic field, also in Fig. [2I^b), show 
that the transition from anti-parallel to parallel orienta- 
tion of the ferromagnetic lead's moments is smooth and, 
moreover, noticeable changes can be observed even for 
fields as tiny as /icxt=0.0001 (a discussion of how small 
this field is in physical units is below). 

The relative orientation of the magnetic moments of 
the leads, and their rotation with magnetic fields, in- 
duce substantial modifications in the conductance and a 
concomitant large magnetoresistance. Figure [31Ja) shows 
the conductance of the trilayer vs. magnetic field ft-ext- 
At fields zero or very small, the anti-parallel orienta- 
tion of the lead's moments produces a very small con- 
ductance. This is reasonable since the spin orientation 
which can conduct in one lead, is blocked by the anti- 
parallel lead. However, as /loxt increases, there are sub- 
stantial changes in the conductance. For fields as small 
as /iext=0.0002, the conductance has changed by about 
3 orders of magnitude already, at least for the particu- 
lar system studied here. Note that if t is assumed to be 
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0.1 eV, then /ioxt=0.0001 is approximately 0.1 T. The 
conductance increases further, by more than an order 
of magnitude, by increasing h^xt toward the value 0.003 
where the moments of the leads become essentially par- 
allel. The conductances remain in absolute value much 
smaller than the perfect conductance 2e^/h (h: Planck 
constant), but its relative changes can be large, as illus- 
trated in Fig. [Sfb). In the widely used definition for con- 
ductance changes, which has the zero field conductance 
in the denominator, the magnetoresistance ratio can be 
very large, and it reaches almost 200,000% at /iext=0.001. 
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FIG. 3: (Color online) (a) Conductance G'(/icxt) in units 
of 2e^/h (h: Planck constant) and (b) magnetoresistance 
(MR) ratio, as a function of an applied magnetic field 
feext. The results are calculated for the ID DE model with 
Jh=8.0, a=1.0, and r=0.0. Positive background charges are 
n^^=n^' =0.65 and n^_^^ = 1.0, and leads and central region 
sizes are L'^'=L'^'=20 and L'°^=4, respectively. To esti- 
mate the MR ratio, defined in the vertical axis of (b), G(/ioxt) 
with /lext = 10"^, instead of G(/ioxt = 0), is used since the 
conductance at /icxt=0 vanishes. Note that with increasing 
/loxt the conductance G'(/ioxt), and the MR ratio, increases 
almost linearly except for very small ftext (^ 0.0003). 



To give a better perspective of the importance of these 
numbers, note that in a recent numerical study2^ of the 
same model used for the trilayer geometry but defined 
on a finite two-dimensional cluster without interfaces, 
the CMR phenomenon was observed for larger magnetic 
fields ft.oxt=0.05. In this case the MR magnitude was 
10,000% at the best. Thus, the trilayer geometry investi- 
gated here certainly produces a more dramatic effect at 
smaller fields than in the bulk simulations. 

To help experimental readers who are more used to re- 
sistance plots, in Fig.m^a), the resistance R of the trilayer 
with an even number of sites in the barrier is shown as a 
function of magnetic field hext- Notice the rapid change 
in R at small /icxt, which is the main result of this paper. 
In experimental similar plots, such as those reported for 
STO and LAO as barriers, hysteresis loops are often ob- 
served in TMR vs. /lext, and moreover the effects at tiny 



fields are negligible .-Si^Jil These differences are caused by 
the anisotropics present in real experiments due to strain 
in the samples, as briefly explained in Sec. II C[ effect 
that has been mainly neglected in our calculations (with 
the exception of the results shown in Fig. [15] below). In 
Fig. H^b) , magnetization M of the classical t2g spins for 
the trilayer ensemble vs. ftcxt is shown. Compatible with 
the spin arrangements already described, there is no net 
magnetization M at /iext=0, since the contributions of 
the ferromagnetic leads cancel out. However, with in- 
creasing /lext, there is a rapid increase in M because the 
lead's magnetic moments are aligning in the same direc- 
tion. The subsequent small increases in Af and decrease 
in R, say for /iext=0.001 or larger, are due to further fine 
alignment of the spins close to the central region (the 
spin orientation in the central region remains antiferro- 
magnetic, which eventually turns to FM with much larger 
/lext)- Thus, the effective height of the barrier melts as 
hcxt (or T, as discussed in the next section) increases in 
this system. 
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FIG. 4: (Color online) (a) Resistance (inverse of conduc- 
tance) vs. magnetic field /icxt for an even number of sites 
^(C)_^ in the central barrier at zero temperature. The leads 
each has I/''^'=I/^'^'=20 sites, and the positive background 
charges are n^|^''=n^|^-'=0.65 and n^^=1.0 in the leads and 
central region, respectively. The couplings used are Jh=8.0, 
and Q=1.0. Notice that the large changes in resistance occur 
at very small magnetic fields, (b) Total magnetization M of 
the classical t2g spins vs. /lext, indicating the rapid develop- 
ment with /lext of a net magnetization. Ms is the maximum 
possible magnetization. 



2. Odd number of sites in the barrier 

The numerical results obtained for the case of an odd 
number of sites in the barrier are very different from those 
reported thus far. In Fig. [51 the results for the case of 
5 sites in the barrier are reported. Fig. [5| (a) 
shows the electronic density n(z) vs. i, which is similar 
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to the results for the case of L(C) = 4 sites in the barrier 
[Fig. [5] (a)]. The electrostatic potential for electrons 
caused by the long-range Coulomb interactions, shown in 
Fig. [5] (b), is also canonical: in the central region more 
electrons are expected to accumulate since there 
is larger than in the leads. The important qualitative 
difference with the previous results in Sec. IIII A H is pre- 
sented in Fig. [5] (c). Here, the spin correlation functions 
for the classical t2g spins indicate that the magnetic mo- 
ments of the leads are parallel to each other, even in the 
absence of magnetic fields. This is in agreement with the 
qualitative scenario for these geometries described in the 
introduction fSec. IIB[) . As a consequence, the case of an 
odd number of sites in the barrier does not present the 
same huge magnetoresistance at small magnetic fields as 
for the even case. 
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FIG. 5: (Color online) (a) Local electron density n{i), (b) 
electrostatic potential </>(«), and (c) classical t2g spin corre- 
lation function (Si • S-,) from the central site at The 
Monte Carlo techniques and optimization method are used 
for the ID DE model with Jh=8.0 and a=2.0 at T=0.01 
and T = 0, respectively. The positive background charges 
are n^-'=n^' =0.625 and n^'' = 1.0, and the system size used 
is L(^'=L'^'=8 and L('^'=5. Here, 4i{i) is self-consistently 
determined. The interface positions are denoted by vertical 
dashed lines. 



The resistance R vs. magnetic field /icxt is shown in 
Fig- Eta). As in the case of the even clearly R de- 

creases with increasing ft-cxt- However, the scales involved 
are very different. While for L^^) even, there are huge 
changes at /loxt as small as 0.0001 (Fig. S]), for the L^^) 
odd case the resistance remains almost the same up to 
/iext=0.02. The reason is that the magnetic field does 
not need to align the magnetic moments of the leads in 
this case (they are already aligned), and moreover the 



up-down-up arrangement of the central region is com- 
patible with the parallel leads' moments. The only mod- 
ification needed in the spin arrangement is the correction 
in the orientation of the central spins that are pointing 
the wrong way. This process takes place between fields 
/ioxt=0.02 and 0.05 approximately. After that, the en- 
tire system is ferromagnetic and the resistance remains 
constant. For completeness. Fig. ^h) shows the total 
magnetization of the classical t2g spins. As expected, 
initially, M is very robust, and it becomes saturated at 
/icxt 0.05. 
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FIG. 6: (Color online) (a) Resistance (inverse of conduc- 
tance) vs. magnetic field /loxt for an odd number of sites 
I/'*^^=5 in the central barrier at T=l/800. The leads have 
L^^^=L''^^ — 8 sites, and the positive background charges are 
n^-'=n^' =0.625 and n^''=1.0 for the leads and central re- 
gion, respectively. The couplings used are Jh=8.0 and a^l.O. 
(b) Total magnetization M of the classical t2g spins vs. mag- 
netic field /loxt, indicating that the system is overall ferro- 
magnetic at all magnetic fields. Ms is the maximum possible 
magnetization. is determined self-consistently at T'=0.2 
and ftcxt = 0, and then used for other temperatures and /loxt 
(see more details in the text). 



B. All-manganite trilayer geometry 
using quantum spins 

The existence of parallel or anti-parallel arrangements 
of the magnetic moments in the left and right leads is also 
investigated for the ID DE model with quantum localized 
t2g spins Si [Eq. dH)] at T = by the DMRG technique. 
For numerical simplicity, spin-1/2 is employed, instead of 
the realistic value 3/2. 

In Fig. [7l the (z-component) spin correlation func- 
tions {SfSj) are reported for a typical set of parameters, 
choosing one of the central-region spins (j = 0) as a ref- 
erence. There are 4 and 5 sites in the central region in 
Fig. [7] (a) and (b) , respectively. In the case of even num- 
ber of central sites [Fig. [7] (a)], at least for large values 
of a such as 2, the antiferromagnetic spin correlations 
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in the central region are clearly observed, and moreover 
the magnetic moments of the left and right leads tend 
to align anti-parallel. However, quantum fluctuations, 
enhanced by the spin-1/2 and one-dimensionality nature 
of the model, make the spin correlations at longer dis- 
tances weaker, and the correlations become very small 
for \i\ ^ 4. Similarly, quantum fluctuations reduce spin 
correlations at large distances also for odd number of cen- 
tral sites, shown in Fig. [7] (b). However, in this case, a 
tendency of parallel alignment of the magnetic moments 
between the left and right leads can still be observed. 
Thus, overall features observed in the case of quantum 
t2g spins are in good qualitative agreement with the ones 
observed in the classical t2g spin case. Realistic two- 
or three-dimensional manganite trilayers are expected to 
behave more similarly to the classical t2g spin case than 
the quantum spin-1/2 one-dimensional case, where the 
quantum fluctuations are the strongest. In fact, quan- 
tum fluctuations appear detrimental for the performance 
of the device proposed in this paper. Higher dimensional 
arrangements will have a stronger tendency to spin order, 
as the classical spins do here. 
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FIG. 7: (Color online) Spin correlation functions (Sf 5|) be- 
tween quantum t2g spins (spin-1/2), calculated from the cen- 
ter of the chain at site j = 0, for the ID DE model with 
Jh = 32 and various a's (indicated in the figure) at T = 0. 
(a): L(^) = i'") = 12 and i^^' = 4, and the positive back- 
ground charge densities are n^'=n^'=0.5 and n^_^^ — 1.0. (b): 
L^^^ = L^^^ = 14 and L'*^^ = 5, and the positive background 
charge densities are n!j^'^=n^'' =0.571 and n^' = 1.0. The nor- 
malization {SoSq) = 1 is adopted. 



With decreasing a, the electron density in the cen- 
tral region is substantially reduced [see Fig. [8] (a)], and 
eventually the antiferromagnetic correlations at short dis- 
tances within the central region disappear for both even 



and odd number of sites in the center (Fig. [7]). Together 
with this effect, an increase in the conductance with de- 
creasing a is observed as shown in Fig. [8] (b). Here 
the conductance is calculated with the time-dependent 
DMRG techniqu o^°i^^ , already explained in Sec. IIIBI 
From these results, in order for the central region to play 
the role of a tunneling barrier with low conductance, it is 
clear that the spins in the central region must be antifer- 
romagnetically aligned, and the electronic density there 
must be close to 1, which is achieved only by a large a, 
i.e., strong Coulomb interactions. 




a 



FIG. 8: (Color online) (a) Electron density in the central 
section ric {— '^n^Qni/L^'^''), and (b) conductance G, as a 
function of a, for various values of the total electron density 
n indicated in the plots. Results are for the ID DE model 
with quantum t2g spin (spin-1/2) and Jh = 32 at T = 0, 
calculated using the DMRG method. The other parameters 
are L'^' = L'-^'> = 12, L'^' = 4 (circles and squares), L'^' = 
= 14, L^*^) = 5 (diamonds), and L*^' = L'^' = 12, 
Z/'*~'' = 5 (triangles). In all cases, the positive background 
charges are n^' = 1.0 and n^^''=n^^\ 



The influence of magnetic fields is also studied in Fig. [SI 
Here, instead of applying an external magnetic field, the 
total z-component of the spins (i.e., total magnetization), 
which is a good quantum number of the system, is var- 
ied as a parameter. Fig. [S] (a) indicates that the anti- 
parallel correlations of the lead's spins, for the case of 
even number L^'-'^ = 4 of central sites, becomes parallel 
with increasing total magnetization. This is in qualita- 
tive agreement with the observation in the previous sub- 
section, indicating that classical and quantum t2g spins 
behave qualitatively similarly in this system. The same 
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FIG. 9: (Color online) Spin correlation functions {StSj) be- 
tween quantum t2g spins (spin-1/2), calculated from the cen- 
ter of the chain at site j = 0, for the ID DE model with Jh = 
32, a=1.0, and various total (i.e., total magnetization) in- 
dicated in the figure. The calculations are performed with the 
DMRG method at T = 0. (a): L^^' = L^^' = 12, L^^^ = 4, 
and total electron density n=0.571. (b): L*^' = L*^' = 12, 
L'*~'-' = 5, and n=0.586. The positive background charge den- 



... (L) (R) 

sities are n . =n\ 



= 0.5 and n 



(C)_ 



= 1.0. The normalization 



(So'S'o) = l is adopted. 



occurs for the case of an odd number of central sites, 
shown in Fig. [5] (b). There is a small caveat to mention 
here: (i) the behavior with increasing total magnetization 
is not monotonous, and (ii) some of the spin correlations 
in Fig. [5] (as well as in Fig. [7]) are small, which might 
be explained by canting effects or by nearly orthogonal 
spin configurations. These issues are not further explored 
here, since they do not appear in the classical t2g spin 
simulations which seem more realistic to describe man- 
ganites. In spite of these caveats, it is clear that quali- 
tatively the similar feature is observed for both models 
with classical and with quantum t2g spins, the key fea- 
ture of having anti-parallel magnetic lead configurations 
for an even number of central sites. 



IV. INFLUENCE OF TEMPERATURE 



up to high temperatures, more specifically room temper- 
ature or above. In fact, previous experimental realiza- 
tions of TMR devices suffered a rapid degradation with 
increasing temperatures,*' as mentioned in the introduc- 
tion (Sec. II A| ). Unfortunately, our proposed trilayer sys- 
tem also has problems in this respect, as shown below, 
but possible avenues to solve this issue are discussed. 

In Fig. [ini classical t2g spin correlation functions for 
two particular sites are shown as a function of tempera- 
ture T, computed with the Monte Carlo technique. One 
of the chosen sites is in the left lead, site 5, and the other 
is in the right lead, site 16, and they are symmetrically 
arranged with respect to the center for the L = 8 -I- 4 -I- 8 
cluster used. For comparison, the same correlation is 
shown for the case when the barrier is removed and the 
entire system is now a unique ferromagnetic metal. The 
latter decays with temperature as shown in Fig. 1101 in- 
dicating that the long-distance ferromagnetic tendencies 
survive up to Tq ~0.03 approximately. Although this 
should not be considered as a critical temperature, due 
to the one dimensionality of the problem which intro- 
duces strong fluctuations at finite temperatures, at least 
it provides a good indicator of the strength of ferromag- 
netism as T is increased. Weak couplings into higher 
dimensional structures will likely stabilize this character- 
istic temperature into a true critical temperatureJ^ii^ 
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FIG. 10: (Color online) Classical t2g spin correlation func- 
tions (S,; ■ Sj) between site i=5 (belonging to the left lead) 
and site j=16 (belonging to the right lead), as a function 
of temperature T. Monte Carlo simulations are used for 



the ID DE model with Jh=8.0, ni^'=n^^> = 



=0.625, nf^ = 1.0, 



l(^)^l(n)^g^ and L'-'^'>=4:. The values of a used are indi- 
cated in the figure. Here, (/!>(i) is determined self-consistently 
at r=0.2 and is subsequently used for other lower tempera- 
tures. For comparison, the classical t2g spin correlation func- 
tions for a uniform DE model with Jh=8.0, L=20, and the 
number of electrons A^=12 are also plotted by circles. 



The large MR effect at low temperature observed in the 
trilayer geometry described in this paper is an interesting 
effect worthy of experimental confirmation. However, for 
practical applications, this large MR effect should survive 



The same spin correlation functions but now in the 
presence of the barrier is also shown in Fig. [TU] for various 
values of a. At low temperatures, these correlation func- 
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tions have the opposite sign (minus) as compared to the 
ones without the barrier. This is simply because in this 
trilayer case there is an even number of sites in the cen- 
tral region, and therefore the magnetic moments between 
the leads align antiferromagnetically fSec. IIIip . Here the 
calculations are carried out by first obtaining the elec- 
trostatic potential at high temperature, where the 
spins are not ordered, and then keeping this potential 
the same as the temperature is reduced. This procedure 
considerably alleviates the numerical effort, particularly 
regarding the Poisson's equation iterations, and tests in 
small systems have shown that this trick does not alter 
the results qualitatively. As also seen in Fig. [TOl the 
spin correlation functions obtained by this method have 
a relatively minor dependence with a. The main point 
to remark is that these spin correlations become negli- 
gible at a temperature T* ^ 0.01, which is considerably 
smaller than the relevant temperature of the pure ferro- 
magnetic system without the barrier. As a consequence, 
it is clear that a strong similarity with previous experi- 
mental results for STO barriers^ may exist in this case. 
Expressed qualitatively, our results indicate that there is 
a temperature scale T*, much smaller than the critical 
temperature Tc, where the orientation of the leads' mag- 
netic moments ceases to be antiferromagnetic. As shown 
below, this is correlated with important changes in the 
magnetoresistance effect. In fact, the large MR effect at 
low temperature seems to occur due to the existence of 
an effective coupling J^s which produces the anti-parallel 
alignment of the leads' magnetic moments. This is rea- 
sonable, since Jeff is an effective weak coupling across a 
tunneling insulating barrier. When the temperature is 
of the order of this coupling or larger, the anti-parallel 
arrangement is no longer preferable and in our model the 
leads' moments rotate freely with respect to one another. 

A systematic study of the influence of the size of the 
central region L'^'^^ on the classical t2g spin correlations 
between the leads is reported in Figure [TlJ As expected 
(Sec. HH), for L(C) even (odd) these spin correlations 
are negative (positive) at low temperatures, indicating 
strong antiferromagnetic (ferromagnetic) correlations be- 
tween the ferromagnetic leads' moments. The magnitude 
of this correlation decreases as L^^^ increases, in agree- 
ment with the expected reduction of the effective cou- 
pling Jos discussed above. Thus, the thinner the barrier 
is, the better the temperature effects become in the tri- 
layer geometry proposed here, namely the higher T* is. 

Our discussion thus far has been based on the spin 
correlations at finite temperatures. It was concluded that 
the large resistance caused by the anti-parallel configura- 
tion of the ferromagnetic leads' moments does not survive 
all the way to the ferromagnetic critical temperatures of 
the leads. This is indeed observed in Fig. [T^ where at 
T* the MR effect is reduced to zero at fields such as 
/Joxt=0. 001-0. 003, while at zero temperature [Fig. [3] (b)] 
the MR effect was huge at the same fields. However, 
even though the MR effect is truly enormous at very 
low temperature, at higher temperatures it is not negligi- 
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FIG. 11: (Color online) Classical t^g spin correlation func- 
tions (Si • Sj) between sites i=5 (belonging to the left lead) 
and site j=12 -I- L^'^-' (belonging to the right lead), for the ID 
DE model with Jh=8.0, q=1.0, n^|^'=n^f' =0.625, n^^°> = 1.0, 
L(l)=L(^)=8, and ^('='=3,4, • ■ ■ ,8. Here, 4>{i) is determined 
self-consistently at r=0.2 and then used for other tempera- 
tures. 



ble, at least in the several Teslas scale of magnetic fields 
/loxt- For example, in Fig. [T^ the MR ratio is plotted 
vs. temperature, for different values of /lext- Note that 
[G(/icxt) — 0(0)] /G(0) is the same as the more standard 
definition [i?(0) — i?(/iext)]/^(^cxt), where the resistance 
i?(/i^^t)=l/G(/icxt)- For fields such as /icxt=0.02 and 0.03 
- corresponding to 20 and 30 T, respectively, if it is as- 
sumed t~l,000 T - the MR ratio can be as large as 100 
or 500% at temperatures of the order of Tc/2. 

In the calculations described so far, no "easy axis" has 
been chosen. In other words, no anisotropics were in- 
troduced. However, the introductory discussion suggests 
that many materials, particularly when in thin-film form, 
do have an easy axis mainly due to the influence of the 
substrate (Sec. II C|) . If due to this effect the magnetic 
moments of the leads cannot rotate "freely" with respect 
to one another as isotropic Heisenberg vectors, but are 
pointing along a particular direction as Ising variables, 
the temperature scale T* above which the magnetic mo- 
ments of the leads become uncorrelated should increase: 
Heisenberg-like isotropic spins can be disordered by ther- 
mal fluctuations more easily than Ising- like spins. On the 
other hand, making more "rigid" the anti-parallel mag- 
netic connection between the leads moments will also pre- 
vent their alignment in very small magnetic fields: larger 
magnetic fields may be needed to achieve the same ef- 
fect as before, i.e., large MR effect. These competing 
tendencies will be discussed in more detail below. 

To investigate this issue, extra couplings between clas- 
sical spins are added to the model, as shown in Fig.[T51 
Via direct Heisenberg couplings between the classical 
spins, Jaf (antiferromagnetic) and Jp (ferromagnetic). 
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FIG. 12: (Color online) Temperature dependence of the MR 
ratio [(G(/icxt) -G(0))/G(0)] x 100 for the ID DE model with 
Jh=8.0, a=1.0, n!j^>=nf =0.625, n^^'=1.0, L'^'^L'^'^S, 
and L^'~'^=4. Applied magnetic fields h^xt are indicated in 
the figure. is determined self-consistently at r=0.2 and 
used for the other temperatures. Note that the effect reported 
at r=0 in Fig. |3] is much larger in magnitude of the order of 
200,000% at /ioxt=0.002 (in that figure the factor 100 was 
not used as in here). In the range of temperatures shown in 
this figure, the original r=0.0 nearly perfect antiferromag- 
netic alignment of the leads' magnetic moments is already 
lost at /iext=0, and thus the changes in resistances are not as 
dramatic as observed at T = 0. 
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FIG. 14: (Color online) Classical t2g spin correlation func- 
tions (Si • Sj) between site i = 5 (belonging to the left lead) 
and site j — 16 (belonging to the right lead), as a function of 
temperature T. Monte Carlo simulations are used for the ID 
DE model with Jh=8.0 and a — 1.0. The positive background 
charge densities are n^:^^ =n^:^'' = 0.625 and nf^' — 1.0, and 
the cluster used is L^^' = L^'^' = 8 and L<*^' = 4. The results 
with Jaf = Jf = 0.1 are shown by triangles, and the results 
with Jaf = 0.1 but no Jf are denoted by diamonds. For 
comparison, the results for a uniform DE model without the 
barrier (circles) and the results with the central barrier but 
no Jaf and Jf (squares) are also plotted. These are repro- 
duced from Fig. 1101 Note that with nonzero Jf and Jaf, the 
spin correlations now vanish at almost the same temperature 
as in the case without the barrier (circles). Here, (/f>(i) is de- 
termined self-consistently at r=0.2 and is subsequently used 
for other lower temperatures. 
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FIG. 13: (Color online) Schematic representation of the 
ID DE model [Eq. 0] for a LSMO/LMO/LSMO trilayer 
(see also Fig. [1]) with additional Heisenberg couplings be- 
tween classical t2g spins, which are introduced to simulate 
anisotropics. Antiferromagnetic Jaf, as well as ferromagnetic 
Jf, couplings are added only in the central region and the 
bonds connecting to the central region, as indicated in the 
figure. The two interface positions are denoted by vertical 
dashed lines. 



both the antiferromagnetic spin arrangement in the cen- 
tral region and the coupling between the center and the 
leads can be made more rigid. That this procedure helps 
regarding the T* problem is clear in Fig. [T3] where the 
classical t2g-spiii correlation functions are shown. Con- 
trary to the case of weakly coupled leads, now a ro- 
bust spin correlation between them survives up to Tq. 
However, the MR effect is not improving at very small 



fields (see Fig. On the contrary, larger values of 

/loxt are needed to achieve the same MR effect as be- 
fore (for instance, the MR ratio at /iext=0.01 without 
anisotropics is 3000% at T-O.OOl, while it is 1500% here 
with anisotropy). Thus, as already mentioned briefly, 
the overall conclusion is that there are two competing 
tendencies to consider in these investigations: (i) on one 
hand, a weak coupling among the leads, mediated by 
the central region, is needed in order for a tiny magnetic 
field to cause a huge effect in transport; (ii) on the other 
hand, for the same reason, a small temperature can en- 
tirely wash out the effect. Adding anisotropics increases 
the effective coupling between the leads' magnetic mo- 
ments, thus helping with the temperature issue (ii), but 
reciprocally, larger magnetic fields are needed to achieve 
the same MR effect. Investigating the subtle balance be- 
tween these competing tendencies is a great challenge, 
both for theorists and experimentalists. 
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FIG. 15: (Color online) The same as in Fig. [12] ex- 
cept that additional couplings between the classical t2g 
spins, Jaf=Jf=0.1, are added to simulate the presence of 
anisotropies (see Fig. I13p . The applied magnetic fields ^ext 
are indicated in the figure. The electrostatic potential (/)(j) is 
determined self-consistently at T=0.2 and then used for other 
lower temperatures. 




FIG. 16: (Color online) Electron density n{i) for the ID 
LSMO/STO/LSMO model (see the text) with Jh = 8.0 and 
a — 1.0, at r = 1/800, and for various values of Vm (in- 
dicated in the figure). The positive background charges are 
(L)_^(R)_g g25 and n!^-'=0.0, and the system size studied 
L(">=8 and L(^'=4. The positions of the two inter- 
faces are denoted by vertical dashed lines. The electrostatic 
potential (/f)(i) is determined self-consistently. 



n 



V. TI-OXIDE BARRIER 

Although the focus of our effort has been on the 
LSMO/LMO/LSMO trilayer, as mentioned in the intro- 
duction we have also carried out model Hamiltonian sim- 
ulations for the more standard case of a non-magnetic 
insulator as a barrier. A widely studied material for 
this purpose is STO, which is here mimicked simply by a 
tight-binding Hamiltonian for the barrier, with its energy 
levels shifted by a one-particle site potential V,„ which 
controls the height of the barrier. In this section, results 
for the case LSMO/STO/LSMO are briefly discussed, 
with the emphasis on qualitative aspects and compar- 
isons with the case of LMO as barrier. The details of 
the model Hamiltonian were already described in Sec.HIl 
and the technique employed is the MC simulation. Once 
again, it should be remarked that subtle effects such as 
"dead layers" are not considered in this study (Sec. lIC|) . 
and their presence may affect quantitatively our conclu- 
sions. 

The local electronic density n{i) is reported in Fig. [TBI 
for the ID LSMO/STO/LSMO model with different val- 
ues of the site potential in the central region. It 
is observed that n(i) in the central region gradually in- 
creases from nearly zero to ^ 0.4 with decreasing I^n, 
indicating that the barrier hight decreases with Vm, and 
consequently the effective coupling between the leads be- 
comes stronger. Note that the lower band of the lead 
(described by the DE model) is centered at — Jh, and 
thus the band of the central region and the lower band 
of the lead are perfectly aligned (i.e., zero barrier height) 
when Vm = — Jh (and a = 0). 



In Fig. [T71 the MR ratio vs. temperature for the 
case Vm — —7.0 is shown. Notice that the MR ef- 
fect becomes appreciable at temperatures comparable 
to the Curie temperature of the individual ferromag- 
netic leads. This is qualitatively similar to the effect 
observed in LSMO/LMO/LSMO for the case where an 
anisotropy was introduced (Fig. [T5)) . which effectively 
caused a stronger effective coupling between the mag- 
netic moments of the left and right leads. Thus, it ap- 
pears that the value Vm — —7.0 chosen for this example 
allows for a robust left-right coupling (as also expected 
from n{i) shown in Fig. I16p . concomitant with the sur- 
vival of a large MR effect up to the Curie temperature. It 
is also observed in Fig. [T7] that because of the strong ef- 
fective coupling between the leads, the MR ratio for small 
magnetic fields (< 0.003) is not as large as in the case of 
the small effective coupling discussed below (Fig. [T5)) . 

Figure [TH| shows results for the same parameters as in 
Fig. [iTl but simply making the height of the barrier Vm 
much larger, i.e., the left and right leads being nearly 
decoupled. Three effects are obvious to the eye: (i) the 
large MR effect appears now at lower temperatures. For 
example, in Fig.[T8]the MR becomes nonzero at T'^0.02, 
while in Fig. [T7] the same occurred at a higher tempera- 
ture 0.03 — 0.04. (ii) On the other hand, increasing the 
height of the barrier much increases the MR ratio at low 
temperatures and small magnetic fields. For example, in 
Fig. [m the MR ratio is about 1,000% at T-O.OOl and 
h^^t = 0.001, while in Fig. [l7] it is about 200% at the 
same T and ft-oxt- This occurs because the barrier is so 
large in Fig. [18] that in the absence of magnetic fields, 
the magnetic moments of the left and right leads are al- 
most decoupled, and thus applying a very small magnetic 
field is enough to align the leads' magnetic moment and 
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LSMO/LMO/LSMO: (1) If the barrier between the mag- 
netic leads is very large (i.e., the effective coupling be- 
tween the leads is small), the magnetic moments of the 
ferromagnetic leads are nearly decoupled (in the absence 
of anisotropics) at hcxt—0, even at low temperatures. 
Thus, the double effect of large barrier and concomitant 
weakly coupled leads causes a large /icxt=0 resistance, 
which enhances the MR effect at low temperatures as 
compared with a lower barrier (i.e., the large effective 
couphng between the leads). (2) However, the large bar- 
rier and weakly coupled leads are fragile upon increasing 
the temperature. 



VI. CONCLUSIONS 



FIG. 17: (Color online) MR ratio [G(/iext) - G(0)]/G(0) x 100 
vs. temperature T for the ID LSMO/STO/LSMO model (see 
the text) with Jh ~ 8.0, a — 1.0, and Vm — —7.0. The posi- 
tive background charges are n^^=n^' =0.625 and n^'=0.0, 
and the system size studied is L'^^'> =L^^'> =8 and L'-'^^=4. Ap- 
plied magnetic fields /lext are indicated in the figure. Here, 
(/!>(i) is determined self-consistently for all temperatures. 



increase conductance, (iii) The dependence of magnetic 
fields on the MR ratio is very mild (compared to Fig. fTT)) . 
namely changing hcxt from 0.001 to 0.03 only alters the 
results by a factor 2 at the most. This is because the al- 
most decoupled leads moments are forced to align in the 
magnetic field direction by a very small hcxt- A further 
increase of /loxt does not change the moments orientation 
at low temperatures. 
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18: (Color online) The same as in Fig. 1171 except that 
-4.0. 



Thus, an interesting and simple picture emerges from 
these qualitative investigations, which are in agree- 
ment with the results shown before for the case of 



In this paper, the use of LSMO/LMO/LSMO as a 
spin-tunnel-junction is proposed. The main difference 
with other previous efforts is in the use of a manganite 
barrier. This would improve the lattice spacing match- 
ing between the constituents, hopefully also alleviating 
the complications found in previous investigations such 
as the infamous interfacial dead layers of STO/LSMO. 
However, the main point of this study is not just the 
all-manganite character of the trilayer, but the antifer- 
romagnetic properties of LMO. For an even number of 
LMO layers, the spin order in the ground state is such 
that the magnetic moments of the ferromagnetic LSMO 
leads are anti-parallel (while for an odd number of layers, 
they are parallel). An anti-parallel-leads configuration 
has a large resistance. But the effective coupling lead- 
ing to this anti-parallel LSMO-moments configuration is 
weak, rendering the ground state fragile. In fact, numer- 
ical simulations show that very small magnetic fields /loxt 
can alter drastically the original ground state at h^xt — 0, 
by aligning the magnetic moments of the ferromagnetic 
LSMO leads. A very large MR effect is observed in this 
transition, at least at low temperatures. Note that in 
arriving to our conclusions a large number of approxima- 
tions have been made, all clearly described in Sec. II CI 
However, we still expect that our theoretical analysis is 
qualitatively correct and may serve as a motivation for a 
real experimental reahzation of the LSMO/LMO/LSMO 
magnetic junction. 

Other effects, such as the infiucnce of temperature, 
were also considered in this study. Together with the 
analysis of STO as barrier, overall trends were iden- 
tified. Very insulating barriers (inducing a very week 
effective coupling between the leads) can lead to low- 
temperature states which are easily destabilized by small 
magnetic fields, causing a large MR effect. However, 
thermal fluctuation rapidly washes out these large effects. 
Reducing the barrier height or introducing anisotropics 
make the original ground state at hcxt = more ro- 
bust with increasing temperature, but this reduces the 
low-temperature MR effect at small magnetic fields. A 
balance between these two tendencies is needed to find 
optimal trilayers for real devices. 
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Note added: While completing the work described in 
this paper, we received two preprints with interesting re- 
lated efforts: (1) Salafranca et ali^ have reported a the- 
oretical study of an all-manganite heterostructure con- 
sisting of FM electrodes and an AF barrier, similar in 
spirit to ours. However, contrary to our proposed sys- 
tem, the chosen barrier^'* is Pr2/3Cai/3Mn03, which has 
CE-type ordering and the same hole doping as the elec- 
trodes, which were chosen to be LSMO with x—l/i. The 
emphasis of Ref. [13 is not on differences between even 
and odd number of central layers as in the present pa- 
per, but on other interesting effects such as the influence 
of the FM electrodes on the spin arrangement of the bar- 
rier. Thus, Ref. H and our efforts nicely complement 
each other. (2) Yu et al?^ have presented experimental 
results for an all-manganite trilayer using LSMO x=0.3 
as electrodes and LSMO a;=0.04 as barrier (the latter 
being almost identical to the LMO barrier theoretically 
proposed in this paper). Those authors report a huge 



TMR ratio of 30,000% at 4.2 K and with bias voltage 
25 mV (our MR ratio is 200,000% at /ioxt=0.001). The 
barrier thickness is 9 atomic layers. Yu et al. assign 
the large MR effect they observed to thermally activated 
magnon resonances inside the barrier. A detail compar- 
ison between theory and experiment will be carried out 
in the near future. 
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